Pairing in spin polarized two-species fermionic mixtures with mass asymmetry 
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We discuss on the pairing mechanism of fermions with mismatch in their fermi momenta due to 
a mass asymmetry. Using a variational ansatz for the ground state we also discuss the BCS -BEC 
crossover of this system. It is shown that the breached pairing solution with a single fermi surface 
is stable in the BEC regime. We also include the temperatures effect on the fermion pairing within 
an approximation that is valid for temperatures much below the critical temperature. 
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I. INTRODUCTION 

In recent times, pairing in degenerate fermi gas of 
atoms [l| has attracted lot of attention. This is the 
outcome of the rapid advancement in the experimental 
techniques to study and manipulate systems of ultracold 
atoms. With these techniques, it is possible to cool and 
trap one or more hyperfine states of an element and con- 
trol the population in each of these states. Furthermore, 
the interaction strength as well as the sign of the interac- 
tion between two components can be tuned over a wide 
range, using the techniques of Feshbach resonance Q. 
When the coupling is weak and attractive, the fermion 
system can be successfully described with the Bardeen- 
Coopcr- Schrieffcr (BCS) theory. The clinching evidences 
are the experimental measurement of the gap energy Q 
and observations of vortices 0]. Typically, in such situ- 
ations the coherence length is much larger than the in- 
ter/particle separation. However, the picture changes as 
the coupling strength is increased. The Cooper pairs are 
more localized and the superfluidity is realized by Bose 
Einstein condensation (BEC) of molecular boson com- 
prising of a pair of fermions. It is further expected that 
such a phenomenon is a cross over between the BCS and 
BEC regimes. These studies have been extended to in- 
clude two fermion species with imbalanced populations. 
In such cases, instead of a crossover, the system is ex- 
pected to show a very interesting and rich phase structure 
with the appearance of exotic superfiuids. These include 
the existence of interior gap superfluidity with one fermi 
surface, breached pairing with two fermi surfaces. It is 
also possible to have inhomogeneous phases like Larkin- 
Ovchhinnikov-Fulde-Ferrel (LOFF) phase wherein the 
Cooper pairs have nonzero net momentum Q , superfluid- 
ity with deformed Fermi surfaces 0] or a phase separated 
state . These exotic phases emerge when pairing occurs 
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between two species whose Fermi surfaces do not match. 
This can happen when the number densities of the two 
species are different or there is a mismatch in their masses 
or both @, i, 0, H0I3. Till date, the two fermion 
species experiments are with the two hyperfine states of 
the same alkali atom forming the condensate [l4[. Very 
recently two fermion species of different masses, lithium 
and potassium, were laser cooled and trapped to degen- 
eracy fl5| . And another recent work reports the observa- 
tions of Feshbach resonances [l(| with the same system. 
Thus achievement of superfluidity with this mass differ- 
ence could be the next frontier of ongoing experiments in 
ultracold fermions. 

In this paper, we attempt to describe such a system 
by constructing a variational ground state explicitly and 
the gap function in the analysis is determined by mini- 
mization of the thermodynamic potential with the con- 
straints of fixed particle number densities for the two 
species. Minimization of the thermodynamic potential 
decides which phase is preferred at the given densities of 
pairing species. This method has earlier been considered 
to describe cold fermionic atoms with equal masses for 
homogeneous (Tc| as well as inhomogeneous pairing [l9l |. 
This method has also been applied to relativistic system 
like cold quark matter and color superconductivity [2f3 |. 
In the present work, we discuss the possible structures 
with homogeneous pairing both with density asymme- 
try as well as mass asymmetry for the two condensing 
species. 

We organize the present work as follows. In section 
ITT1 we discuss the ansatz for the ground state and the 
Hamiltonian in terms of the Four fermi point interac- 
tion to model the superfluidity for the two species of 
fermionic atoms. In section IIIII we evaluate the ther- 
modynamic potential by minimizing the thermodynamic 
potential with respect to the functions in the ansatz for 
the "ground state" . In section IIVI we discuss the results 
regarding asymmetric fermionic populations and the gap- 
less phases. Finally we summarize and conclude our re- 
sults in section W\ . 
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II. ANSATZ FOR THE GROUND STATE AND 
THE HAMILTONIAN 

To examine the superfluidity for fermionic atoms, 
we consider a Hamiltonian describing two interacting 
fermionic species with four-fermion point interaction 
given as 



h 2 v 2 

2rrij 



**(z) 



(1) 



III. EVALUATION OF THERMODYNAMIC 
POTENTIAL AND THE GAP EQUATION 

Having defined the ground state as in eq.((3]), we next 
evaluate the thermodynamic potential corresponding to 
the Hamiltonian given in Eq.(p}. To calculate, e.g., the 
energy density, one can take the expectation value of the 
Hamiltonian. Noting that the variational state in Eq.([3]) 
arises from successive Bogoliubov transformations, one 
can calculate the expectation values of the various oper- 
ators. Thus we have, with (O) representing the expecta- 
tion value of an operator O in the new ground state of 
the system (Q, (3, fi)\0\Q, (3, fi), 



where r and s are the spin indices and i denotes the 
species with mass m^. The constant g is the bare inter- 
action strength between the two species and is related 
to the s-wave scattering length a. To describe pairing 
between two different fermionic species , we consider the 
ansatz for the ground state [H, Qj| of the system as 



o(B r -B) 



|0>, 



(2) 



where B f = / dke y ¥* t (k)/(fc)*- r (-k)- Here e ij is the 
Levi-Cevita tensor, with i and j denoting two differ- 
ent fermionic species. The function f(k) is the varia- 
tional function related to the order parameter, as will 
be seen later. In the case of equal population and for 
negative weak coupling, this ansatz corresponds to the 
standard BCS wave function. The two ground states, |0) 
and |fi), are related by the unitary transformation oper- 
ator U = e^ B ~ B \ Hence the field operators transform 
as \&' = U^W where as is the annihilation operator 
for |f2). To include the effect of temperature and den- 
sity, we use the method of thermo-field dynamics (TFD) 
that is particularly useful while dealing with operators 
and states. Here, the thermal "ground state" is obtained 
from O) through a Bogoliubov transformation in an ex- 
tended Hilbert space associated with thermal doubling of 
operators. Explicitly, \Q,(3,fI), the ground state at finite 
temperature and density is given as 



where, 



|n,AM>=exp(Sj >p -B / ,, M )|n> 



Bl M = J [*'(k)t0L(k,/?,M)*a-k)] dk. 



(3) 



(4) 



In Eq.Q, the function as we shall see later, will be 
related to the distribution function of the i th species and 
the underlined operators are the operators in the ex- 
tended Hilbert space associated with thermal doubling. 
All the functions in the ansatz in Eq.([3]), the condensate 
function /(k), the thermal functions l (k, (3, fi) shall be 
determined by extremizing the thermodynamic potential. 
We carry out this extremization in the next section. 



(¥jt(ki)*;(ka)) = [cos 2 (/(A :i ))sin 2 (0 1 (fc 1 ))+sin 2 (/(fc 1 )) 
cos a (0 2 (*i))]a„<y(ki-k a ) (5) 

<* 2t (k!)vl/ 2 (k 2 )) = [cos 2 (/(fc 1 ))sin 2 (0 2 (fc 1 ))-sin 2 (/(fc 1 )) 
cos 2 (0 1 (fc 1 ))]<5 rs( 5(fe 1 -fc 2 ), (6) 

<^(k!)*2(k2)) = 



, 2/1 . N \ . sin(2/(fcx)) M ^ 2 



[l-sin^^)) 

- am 2 (e 2 (k 1 ))]5 r ^ s 8(k 1 + k 2 ), 
<^t (kl )^ t(k2)) = sin(2/( fel )) [l _ shi 2 {6l{kl)) 



Note that the thermodynamic potential is given by 
fi = e - -p ~ ViP , 



(9) 



where, e = (H)p tli = T + V is the energy density, with T 
and V as the kinetic and potential energy contributions 
respectively, s is the entropy density and /z, is the chem- 
ical potential for the species 'i'. The diagonal part of the 
potential in Eq. ([9]) is given by 



^3 J d 3 fc[^[cos 2 (/)sin 2 (e 1 ) 



(27T) 

+ sin 2 (/)cos 2 (0 2 )] +e 2 [cos 2 (/)sin 2 (0 2 ) 
+ sin 2 (f) cos 2 (6,)}], (10) 

where, ^ = h 2 k 2 /2rrii — /Ltj is the kinetic energy with 
respect to the chemical potential, of the i th species. Sim- 
ilarly the expectation of the term H lnt is simplified using 
the Wick's theorem 

V=(H int )=gp lP2 +gI 2 D , (11) 

Id is related to the condensate defined as Id = 
5 r - s (k)^ (-k)) and is given as 



1 



(2tt) 3 



dk sin 2/(k) (cos 2 0i - sin 2 6 2 ) . (12) 



(7) 



sin 2 (0 2 (fc 1 ))]5 r _ s <5(fe 1 + fe 2 ). (8) 
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Further, the species densities pi = ^^(k)^* (k)) 5 rs for 
the fcrmions are given as 



P2 



d 3 k 



1 



Pi 



f>2 



d 3 k [cos 2 (/)sin 2 (6»i 



(27T) 3 

+ sin 2 (/)cos 2 (# 2 )] 
1 



(27T) 



l + -)sin 2 (0 2 ) 

UJ 



(19) 



d 3 k \cos\f)s\u z {e 2 ) 



(2tt) 3 
+ sin 2 (/)cos 2 (0 1 )] 



(13) 



(14) 



where, we have denoted £ = e — v. The other varia- 
tional extremisation conditions 5(Q)/S6i(k) = deter- 
mine Oi{k) to be related to the distribution functions for 
the ith fermion species as 



Finally the entropy density for the two species fermionic 
mixture is 1211 



sin 2 (6»i(fc,/i)) 



1 



exp(/?cj. 



1 



(20) 



1 



d 3 k[n l (k)\og{n i (k))+ 



^ (2tt) 3 
(l-n i (fc))log(l-n i (fe))] J 



(15) 



In the above, the quasiparticle energies are given as u>i = 
lo + S^ and uj 2 — uj — 5% where 5% = [(ei — v\) — (s 2 — v 2 )]/2. 

Next we examine the gap equation Eq. (Tl"T)) . which for 
nonzero 8 can be written as 



1 1 f 1 

where rii(k) = sin 2 (#i) is the density distribution of the = ,„ ^ 3 / d 3 k— [l — sin 



^O-sin 2 ^) 



(21) 



species. Combining Eqs.tflO]), (JTTJ, and one can 
then calculate the expectation value of the thermody- 
namic potential in the ansatz state given in Eq.([3|. The 
thermodynamic potential is a functional of three func- 
tions, the condensate function, f(k) and the two thermal 
distribution functions 6i(k) for the two species. These 
functions arc determined by functional minimization of 
the thermodynamic potential of Eq.@ which we shall 
analyse in the next subsection. 



A. Gap equation 

Functional minimization of the thermodynamic poten- 
tial f2 with respect to f(k) gives 



tan(2/(fc)) = 



(ei + e 2 ) - iy x + v 2 ) 



(16) 



where z/j = pi—ge tJ \pj\ is the chemical potential with the 
mean field correction. We have also defined in the above, 
the superconducting gap A = — gljj and e = (ei + £2)/2, 
v = (y\ — v%)/2 as the average kinetic energy and chem- 
ical potential respectively. Let us note that, the con- 
densate function depends on the average kinetic energy 
and the average chemical potentials of the two condens- 
ing species. Substituting the solution for the condensate 
function from Eq. (Tl6|) in the definition of Id given by 
Eq. (fT2"j) . we obtain the gap equation 

A = ~ j£r I d3h ^ [1 ~ sin2(0l) ~ sin2(02)] ' (17) 

Similarly, one obtains the densities for the two fermion 
species as 



1 



Pi 



(2tt) 3 



d 3 k 



i(l-i)cos> 2 



(18) 



This equation is ultraviolet divergent which is character- 
istic of the contact interaction. It is rectified by subtract- 
ing the contribution at T = and /i = and relating this 
renormalized coupling to the s-wave scattering length a 
[13, EH • Thus regularized gap equation is 



Anh 2 a (2?r) 



1 ffkfai-rtw 



sin 2 (# 2 )] - — ) , 



(22) 



where rh is the reduced mass. 



B. Stability condition 

The stability of the pairing state is decided by compar- 
ing the thermodynamic potential of the superconducting 
matter with that of the normal matter. Thus the relevant 
quantity is the difference of the thermodynamic potential 
between the paired and normal phases. The thermody- 
namic potential of the paired fermionic mixture is 



O 



(27 



If 1 

— J d 3 k £-w--^ln(l + (*p (-#*)) 



A 2 
9 



9P\Pi- 



(23) 



Subtracting the thermodynamic potential for normal 
matter (A = 0) from the above equation, we have the 
difference in the thermodynamic potential between the 
condensed and the normal matter as 



sn 



(2 



^3 J d 3 k |e| - w - i^ln (1 + cxp (-#*)) 

(24) 



--^ln(l + exp(-/?w i)) 



A 2 
9 
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Here a;, = u ± 5^ and ujoi = |£| ± <5^. This difference in 
the thermodynamic potential, SQ has to be negative for 
the stability of the paired state. Further one can use the 
gap equation to eliminate the coupling g in Eq. (|24[) to 
obtain 



so 



(2tt) 3 



i 

exp (-/Jwi)) + + ex P (-/ ?CJ 0i:)) 



2w 



[sin 2 (00 



sin' 



(25) 



This expression is free of any ultraviolet divergence and 
will be used to determine the stability of the given paired 
state. 



£1 



x 2 dx 





4) 







(l-0(-£ 2 )) 



(31) 



x dx 





■1) 







(i-e(-^)) 



(32) 



The Eqs. ([30l) - (|32|) are the governing equations of the 
paired state for two component fermionic mixture at zero 
temperature, when the components arc homogeneous and 
not equal in density. These three equations are to be 
solved self consistently. In the zero temperature limit, 
the difference in the thermodynamic potential 



C. Zero temperature limit 

In the limit of zero temperature, the quasi-particle den- 
sity distribution of the two atomic species 



rii(k) = lim 



1 



/3^oo exp(/3wi) + 1 



0(-w<), (26) 



where 0(. . .) is the Heaviside step function. The gap 
equation in this limit is given as 



1 



(2tt) 3 
1 

~2Tk 



d A k 



1 



_(i-e(-wi)-e(-wa)) 



(27) 



The densities of the two atomic species are 
1 



Pi 



P2 



(2tt) 3 



d 6 k 



i 



(2tt) 3 
1 



d d k 



1 



(l-e(-wi)) 



(28) 



(29) 



For numerical calculations, it is useful to express the 
equations Eqs. |25|) and ([27|) ~([29 |) in terms of dimcn- 
sionless quantities. Hence we make the substitutions k = 
Ufx, A = £j?A, v = epp, 8 V = ep8 w and u> = tpCj where 
kp is the Fermi momentum defined as kp = 3ir 2 (pi + p 2 ) 
and ep = h 2 kp/2m. In terms of the dimensionless quan- 
tities 



7T 

kpa 



i - e(-wi) - e(-wa) - \ 



(30) 



d A k 




A 2 



1 



(27T) = 



d A k 



Q(—CJl) + ( U>2 



d A k 



A 2 \ 

5j) e ( -^) 



W l6(-W0l) +W O 20(-Wo 2 ) 



(33) 



In the above, u>\ = u) + 5^, L02 — u> — 6^, loqi = |£| + 5% 
and uj 02 = |£| - 5(. 



D. Breached Pair solution 

For the symmetric case of two fermionic species of 
equal masses, having equal densities, 6% is zero. Then, 
= W2 = uj and the Q(—u>i) in the equations are zero 
as uji > always. This is the standard BCS phase. In the 
general case when 6^ is nonzero, one of the 0(. . .) func- 
tions in the equations has a nonzero contribution. This 
can occur when there is a difference between the densities 
or a difference in the masses of the two fermion species. 
Without loss of generality, let us consider the case when 
6^ is negative. In this case a->2 = ui — is always pos- 
itive and (—LU2) will be zero. However, cji(= uj + 5%) 
is negative when to < \S^\ and in this domain (— oji) is 
nonzero. The quasi particle excitations become gapless 
at momenta fe max and fc m i n given as 



2;„2 

max/min 



h 2 k 



{m\V\ -I- m 2 ^2) 



± \J (m\V\ — miV2) 2 — 4mim2A 2 . (34) 



The fermionic mixture supports a gapless mode in the 
momentum domain (^m!^max)i m scaled units x m i n = 
^min/^F and x max = k max /kp. 

It is convenient to reduce the density equations to the 
total and the difference of the densities. Using these, one 
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can define the experimentally relevant parameter polar- 
ization P = (p± — p2)j{p\ + P2), the measure of popula- 
tion difference of the two fermionic species. Since only 
0(— lo\) is nonzero, the difference of the densities 



27T 2 ( Pl - P2 ) 



hp 



x dx 



;(4 



x min)- (35) 



Rearranging the above equation, we get 

1.3 

s ^ p 1 3 3 \ 

VP — Pi ~ P2 — ~Q^2\ X max ~ 'Z'miiJ- 

The total density in dimcnsionlcss units 

4 = ^/^[e(-^) + (i- 

(l-e(-wi)) 



(36) 



(x 2 /2 - i>) 



(x 2 /2-£>) 2 + A 2 



(37) 



To solve Eqs.([30 |) - ([37|) . the parameters of the interacting 
fermions: scattering length a, density of atoms p and 
masses of the atoms m 2 ; are to be chosen appropriately. 
This can be achieved using the dimcnsionlcss quantity 
kpa with the definition p/kp = 1/3tt 2 . Then the total 
density 



x dx 1 



x dx 



x 2 /2 - v 



(x 2 /2-v) 2 + A 2 
x 2 /2-P 



(x 2 /2-v) 2 + A 2 



(38) 



For the value of p/kp defined earlier, the polarization 



P 



(39) 



This is an important relation as P is a control param- 
eter in experiments, which can be varied over a wide 
range with fine control. Further, to study the mix- 
ture of fermionic atoms with unequal masses, we de- 
fine the dimensionless quantities a = (m\ + 7712) /2rh, 
j3 = (mi — m2)/2m, and 7 = mim2/m 2 . In terms of 



these quantities k 
as 



max/mm 



= k 



F x v 



max /mm 7 



which arc given 



max/ mm 



(av + P8 V ) ±J(/3P + a8 u ) 2 - 7A 2 . (40) 



While studying mixtures of different mass ratios, it is 
also convenient to use the mass ratio q = mi/7712 as a 
generic parameter. To study finite temperature effects, 
temperature is expressed in units of fermi temperature 
Tp, defined as ksTp = ep. 



(kpa) 



FIG. 1: A is plotted against coupling (kpa) 1 at zero tem- 
perature without population imbalance. 



IV. RESULTS AND DISCUSSIONS 

Let us note that the dimensionless parameters which 
describe the fermionic mixture are the dimensionless cou- 
pling (fcpa) -1 , polarization P, temperature T in units of 
Tp and mass ratio q. The gap equation, Eq. (|30p . to- 
gether with the number density equations, that is, the 
average density given by Eq. ([38|) and polarization given 
by Eg. (|3"9"]) , are solved self consistently to obtain order 
parameter A, and chemical potentials v and 5 V . The 
thermodynamic potential, Eq. (j2"4")l is then calculated 
using these parameters. 

We first consider the zero polarization case. Figure [T] 
shows the variation of A with coupling (kpa)" 1 . In the 
BCS regime where (kpa)" 1 < —1, A is exponentially 
small in agreement with analytic solution as given in ref. 
[TgL [22j . The variation of chemical potential v is shown 
in the figure [2l The chemical potential v decreases as the 
coupling (kpa) -1 increases and it is zero at (kpa)" 1 w 
0.55. It is negative at higher values of (kpa)" 1 , which 
indicates the formation of Bosc-Einstcin condensation of 
diatomic molecules of fermionic atoms. 

Next we consider the case of nonzero polarization, how- 
ever, for symmetric mass case i.e. with mi = TO2. In such 
a case, the breached pair phase shall have fermionic gap- 
less modes. The gapless modes occur when to* = 0. With- 
out loss of generality, we shall assume here 5^ as negative. 
For mass symmetric case, q = 1 and hence <5j = —8 V . In 
this case only u/± can become zero. Thus gapless mode 
means lu = \5^\. Depending on the chemical potentials, 
breached pair solution can exist cither with one (yi < 0) 
or two (2^1, v 2 > 0) fcrmi surfaces referred to as BP1 and 
BP2 phase respectively Q. 

The figure [3] shows the density profile for mass sym- 
metric case at zero and finite temperatures in the BEC 
regime of interaction. The density profile here corre- 
sponds to (kpa)" 1 = 2 and P = 0.2 and has the char- 
acteristic of BP1 phase of a single fcrmi surface. We 
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(k F a) 



FIG. 2: v is plotted against coupling (kFa) 1 at zero temper- 
ature without population imbalance. 
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FIG. 4: The momentum space density profile for mass sym- 
metric case q = 1 at (a) T = (b) T = 0.005T F . Here 
(kpay 1 = and P = 0.2. The profile corresponds to BP2 
phase. This phase, however, is unstable. 



mass ratio q = m\jmi differing from unity. Specifically, 
we have taken the mass ratio q = 0.15. This ratio cor- 
responds to 6 Li- 40 K mixture, which has been cooled to 
degeneracy recently 15|, with 6 Li chosen as the majority 
population. The momentum space density profile for this 
system at (hp a) -1 = 0.1 and P = 0.2. is shown in the 
Fig. O Though the coupling strength is close to unitar- 
ity, the phase corresponds to the gapless modes with one 
fermi surface (BP1). The critical polarization turns out 
to be P c w 0.35. 



FIG. 3: The momentum space density profile for mass sym- 
metric case q = 1 at (a) T = (b) T = 0.005T F . Here 
(kFa)^ 1 — 2 and P — 0.2. The profile corresponds to BP1 
phase. 



also verify here that the thermodynamic potential differ- 
ence between the paired phase and the normal matter 
is negative indicating its stability. The critical polariza- 
tion P c up to which this phase is stable for this cou- 
pling is P c — 0.52. As the coupling increases, the critical 
polarization P c increases and finally reaches P c = 1 at 
(kfa)^ 1 « 2.3. The upper and lower curves correspond 
to zero temperature and T/Tf = 0.005. Finite tempera- 
ture effects smoothen the distribution functions. 

We also observe that, at unitarity ((fc/a) _1 =0), 
breached pair solution exists with two fermi surfaces 
(BP2). However it is thcrmodynamically unstable i.e. 
S£l is positive. The density profile in momentum space 
is shown in the Fig. 2J We have also shown the effect 
of temperature in the density profiles within the present 
mean field calculation. As before, the density profiles 
get smoothened for finite temperatures as quasi-particle 
density distribution is no longer a 0(. . .) function. 

We next consider the mass asymmetric case with the 
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FIG. 5: The variation of A against polarization is shown 
for mi/iri2 = 0.15 at couplings (fcpa) -1 = 0.1 (dashed) and 
(kpa) -1 = 0.5 (dot-dashed). 



We see the breached pairing solution here with one 
fermi surface in the deep BEC regime. For example for 
(fcpa) -1 = 8, we find the stable BP1 state up to the 
critical polarization P c ps 0.22. 
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FIG. 6: The momentum space density profile for mass asym- 
metric case q = 0.15 at (a) T = (b) T = 0.005T F . Here 
q = mi/rri2, (/cf^)" 1 = 0.1 and P — 0.2. The profile corre- 
sponds to BP2 phase. This is stable phase. 

V. SUMMARY AND CONCLUSIONS 

We have considered here a variational ground state for 
the system of nonrelativistic fermions with a four fermion 
point interaction to model the phase structure of the 
cold atomic mixture near the Feshbach resonance. The 
ansatz functions including the thermal distribution func- 
tions describing the variational ground state are deter- 
mined through the minimization of the thermodynamic 
potential. The stability of the solutions is decided by 
comparing the thermodynamic potentials of the paired 
state and normal matter. 

We find that gapless modes with a single fermi surface 
solutions are possible within in the BEC region of the 
coupling constant when there is a difference of number 
densities of the two species. We have in particular looked 



into the number density difference of the two species dif- 
fering in their masses. There is a critical polarization for 
a given coupling which can sustain pairing. Beyond this 
value for polarization the pairing state becomes unstable. 
We have not calculated the Meissner masses [23| or the 
number susceptibility [24| to discuss the stability of dif- 
ferent phases. Instead, we have solved the gap equations 
and the number density equations self consistently and 
have compared the thermodynamic potentials. In certain 
regions of couplings, we have multiple solutions of the gap 
equations. In such cases we have chosen the one which 
has the least value for the thermodynamic potential. 

The present results however are limited by the simpli- 
fied ansatz as considered here in terms of fermion conden- 
sates. To look into the thermal effects we have consid- 
ered rather small temperatures. At higher temperature, 
particularly near the critical temperature, the effects of 
fluctuations involving corrections arising from collective 
modes will play an important role, particularly for strong 
coupling. Further, considering a realistic potential rather 
than the point interaction as considered here as well as 
the calculation of some of the transport properties in the 
different phases of cold atomic gas, will be very interest- 
ing to investigate and can be studied within the present 
framework. Some of these calculations are in progress 
and will be reported elsewhere. 
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